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We examine ion acceleration by irradiating a hundred TW laser pulse on a double- layer target. It 
is shown analytically and by three-dimensional particle-in-cell simulations that higher energy protons 
are obtained by using material with a high charge-to-mass ratio in the first layer of a double-layer 
target, because a strong Coulomb explosion occurs in such a material. As a result, the protons keep 
accelerating for a longer time. Using the optimal conditions for the target, it is shown that high 
energy and high quality protons can be generated. 
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Laser driven charged particle acceleration is one of the important examples of applying compact laser technology. 
This method of fast particle generation is very attractive, since the acceleration rate is markedly higher and the 
facility size can be substantially smaller than that of conventional accelerators. Laser driven fast ions are promising 
in many applications such as hadron therapy [1^, fast ignition of thermonuclear fusion [2 , production of positron 
emission tomography (PET) sources [3, conversion of radioactive waste [4 , proton imaging of ultrafast processes in 
laser plasmas [5 , a laser-driven heavy ion collider [6] , and a proton dump facility for neutrino oscillation studies [7l . 

Although there are many experimental studies of laser acceleration of ions (see review articles [8 ), the achieved 
proton energy at present is not high enough for some applications such as hadron therapy which requires one-two 
hundred MeV protons. Higher energy protons can be obtained by using a higher power laser. However, the laser power 
enhancement will result in the cost increase of the accelerator. Therefore, it is important to study the conditions for 
generating higher energy protons with a lower power laser by using some special techniques. For example, one can 
obtain the required energy protons by using radiation pressure dominant acceleration (RPDA) and its modifications 
[61I9UTT]. The laser peak intensity is /q > lO^^W/cm^ in those studies. 

In this Letter, we show a method to obtain the high energy protons by using a laser whose intensity is 
Jo ~ lO^^W/cm^ and energy is £ias ^ 20 J, when the RPDA regime is not reached in full scale. We use three- 
dimensional (3D) particle-in-cell (PIC) simulations in order to investigate how high-energy and high-quality protons 
can be generated by a few hundred TW laser. We study the proton acceleration during the interaction of a laser pulse 
with a double-layer target consisting of some high-Z atom layer and a hydrogen layer (see Fig. [ija)). As suggested 
in Refs. [U [12], a quasimonoenergetic ion beam can be obtained using targets of this type. Our aim is to obtain 
a high energy {£ > lOOMeV) and high quality (Af /f < 5%) proton beam using a relatively moderate power laser 
(P ^ 500TW). We show the dependence of the proton energy on the material of the first layer. In experiments of 
laser driven ion acceleration, the case of using a CH polymer target observed higher energy protons [13 than the case 
of using a metallic target [14 . 

Here we estimate the required laser energy to generate protons with energy £ > lOOMeV. The electric field of a 
charged infinite sheet is equal to Eq = p//2eo, where p is the charge density, / is the sheet thickness and eo is the 
vacuum permittivity. The x component of the electric field of the positively charged disk is equal to 



where R is the charged disk radius. We assume that the x axis is normal to the disk surface placed at the disk center. 
The energy gain of protons in this electric field is £p = qeEx{x)dx = qeE^R^ where Qe is the electron charge. 

Let Ne electrons be extracted from the initially neutral charge disk. The resulting charge density is given by 
p = QeNe/TrR'^l, where / is the disk thickness. The energy of the proton is given by £p = q^N(>/2eQT:R. We obtain 

7V, = ?^^. (2) 

This expression shows the number of electrons that should be removed from the disk to generate an electrostatic field 
which accelerates protons up to the energy of £p. Next, we estimate the energy necessary to remove electrons from 
the initially neutral charge disk provided that the remaining ions do not move significantly. Let N electrons be already 
removed from the disk. The charge density of this disk is, p{N) = qeN/nR'^l, and Eo{N) = p{N)l/2eo = q^N / 2ei^7i R^ . 
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FIG. 1: (a) The first layer and proton bunch shape of the disk case, (b) The initial shape of the target and the exploded first 
layer by the Coulomb explosion, (c) The electric field component, Ex, normalized by Eq, of the disk case (blue solid line) and 
the Coulomb explosion case (red solid line), (d) The energy gain of protons, 8p, normalized by the maximum of the Coulomb 
explosion case, the disk case (blue dash line) and the Coulomb explosion case (red dash line). 



The energy necessary to remove one electron from the disk is Ee{N) = qeEQ{N)R = q'^N/2eo7rR. If electrons are 
removed by the laser pulse, the required energy of the laser pulse is 

£las= / £e{N)dN= (3) 

Jo Qe 

This expression shows the minimum laser pulse energy needed to generate a proton with the enegy of £p in the 
electrostatic field of the disk. When we assume the diameter of disk to be 10/im, the necessary laser energy to 
generate the proton of fp=200MeV is calculated to be Sias ^6 J. This estimation assumes that the energy of the laser 
is efficiently spent on extracting electrons while ions remain at rest. 

In the model described above, the protons are accelerated in the electric field given by Eq. ([T]). This electric field 
decreases rapidly with the distance from the disk surface as shown in Fig. [TJc) (blue solid line). To enhance the ion 
acceleration, we must produce the condition where the ions are being accelerated for a longer time by the electric field 
which more slowly decreases with distance as seen by the protons. This can be achieved using a moving electrostatic 
potential, associated with the motion of the first layer of the double-layer target. With a sufficiently high power laser 
the first layer can be pushed forward in the RPDA regime However, as suggested in a Coulomb explosion 

of the first layer can produce a moving electrostatic potential, by which the second layer is accelerated. 

Here we estimate the effect of this Coulomb explosion of the first layer. To simplify the discussion, we assumed 
that the disk target expands in the shape of cylinder. It is assumed that the target thickness, of the initial shape 
expanded to -\-x and —x and the target radius, expanded to r{x)R^ where r{x) > 1, and the distance of a proton 
from the expanding target surface is x' . The electric field Ex at the proton position Xp — X I X IS 




where Eq is the value for the initial disk. When the r{x) is a linear function of x of the form r{x) = ax + 1, 
a is a real number and a > 0, and the proton is on the expanding target surface, this electric field is Ex{x) = 
Eo{l + aR/2 + R/2x - y^l + (ai?/2)2 + aR'^ /2x + {R/2x)'^}/{ax + 1)^. Here E^ix) and £p (x) increase with the 
decrease of a. In the limit of a = 0, r{x) = 1 (Fig. Bb)), this electric field is E^ix) = Eo{l + R/2x - y^l + {R/2x)'^). 
This electric field Ex{x) is shown in Fig. [TJc) (red solid line). The electric field of this case is higher than the disk 
case at all positions. The energy gain of protons is Sp{x) = qeEo Jj^^^ Ex{x)dx. The proton energy of the Coulomb 

explosion case is 1.6 times higher than the disk case until the position x/R=S (see Fig. [ijd)). In a spherical expansion, 
on the other hand, the proton energy is the same in the expanding target and non-expanding one. 

The acceleration rate is higher when the front of the expanded first layer has the velocity close to the accelerating 
protons. That is, we should make a strong Coulomb explosion of the first layer. The higher the Qi/rni of the ions, 
the greater their energy. A stronger Coulomb explosion occurs in the high qi/rui material. Therefore, we can obtain 
higher energy protons by using the high Qi/rrii material in the first layer. This consideration is corroborated by the 
simulations described below. 
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FIG. 2: (a)-(e) Particle distribution and electric field magnitude (isosurface for value a = 2); half of the box of the electric 
field is removed to reveal the internal structure, (a) Initial shape of the target and the laser pulse, (b),(d) the interaction of 
the target and laser pulse, and (c),(e) the first layer shape and the accelerated protons (color scale). (b),(c) The gold case and 
the (d),(e) carbon case, (f) Distribution of gold ions and protons of the gold case, and (g) the carbon case; a 2D projection is 
shown looking along the z axis. 



The simulations were performed with a 3D massively parallel electromagnetic code, based on the PIC method 
[15]. The dimensional quantities are given in terms of the laser wavelength A = 0.8/im; the spatial coordinates are 
normalized by A and the time is measured in terms of the laser period, 27r/a;. We use an idealized model, in which a 
Gaussian p-polarized laser pulse is incident on a double-layer target. The laser pulse with the dimensionless amplitude 
a = qeEo/rrieOJC = 50, which corresponds to the laser peak intensity of 5 x lO^^W/cm^, is lOA long in the propagation 
direction and is focused to a spot with size 4A (FWHM), which corresponds to the laser power of 620TW and laser 
energy of 18J. Both layers of the double-layer target are shaped as disk. The first layer has the diameter of 8A and 
thickness of 0.5A. The second, hydrogen, layer is narrower and thinner; its diameter is 4A and thickness is 0.03A. 
Although, the first layer material is varied, it is assumed that it is comprised of ions with the charge state of Zi = -\-6 
and the number of ions is the same in all cases. The electron density inside the first layer is ng = 3 x lO^^cm"^ and 
inside the hydrogen layer is ng = 9 x lO^^cm"^. The total number of quasi-particles is 8 x 10^. 

The laser pulse is normally or obliquely incident on the target. The oblique incidence of the laser pulse is realized by 
rotating the target about the Z axis. In the normally incident cases, the number of grid cells is equal to 3300 x 1024 x 
1024 along the X, F, and Z axes. Correspondingly, the simulation box size is 120A x 36. 5A x 36. 5A. In the obliquely 
incident cases, the dimension of Ly is set to be two times larger than the normally incident cases [16]. The boundary 
conditions for the particles and for the fields are periodic in the transverse {Y^Z) directions and absorbing at the 
boundaries of the computation box along the X axis. In figures [2l3l5l the X axis denotes the direction perpendicular 



to the target surface and the y axis is parallel to the target surface, while the direction of z axis is the same as the 
direction of the Z axis. The origin in the coordinates is located at the center of the rear surface of the initial 
first layer. 

We show two results, one is a case using gold for the first layer and the other case uses carbon. In order to examine 
the dependence of the proton energy £p versus the first layer material, we performed simulations with the normal 
incidence of laser pulse. Figures [2|b),(c) show the case of gold for the first layer and Figs. [2|d),(e) show the case of 
carbon. We see that the carbon ions are distributed over a wider area than in the gold case by the Coulomb explosion. 
Figures [2|f),(g) show a 2D projection onto the {x^y) plane. Here we present the region of z = —0.7 ~ 0.7A to see the 
ion density around the z = plane. 

In the carbon case, the center of the exploded first layer moves in the laser propagation direction (Fig. |2|g)). This 
velocity is about 1/10 of that of the protons. We show the estimation of this moving first layer effect on the proton 
energy. We assume that the first layer velocity is V and the second layer ion, proton, velocity is v. That is, the 
accelerating electric field moves at speed V, it is denoted Ex{x — Vt). For mV'^ <C Spo^ we obtain the proton energy 



£p = £pO + V^/2m£po, (5) 



where £po is the proton energy in the case of a non- moving first layer, and m is the proton mass. For the proton 
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FIG. 3: (a) The accelerating electric field, Ex, that the protons receive and the energy gain of protons 8p. (b) The proton 
energy spectrum of the gold case and the carbon case at t=100, normalized by the gold case maximum peak. (c)(d) The spatial 
distribution of particles and the electric field component Ex (color scale) and the electric field component Ex on the x axis at 
t = 50. 
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FIG. 4: Proton energy normalized by the energy of the gold case vs the first layer material. 



velocity in the non-moving first layer vq <^ ewe obtain £p = Spo{l + 2V/vo). V/v ~ 0.1 in our simulation, v/vq ~ 1.1, 
so we obtain 20% higher energy protons by this moving first layer. 

Figure [sja) shows the simulation results of the accelerating electric field that the protons see in the gold and carbon 
case, and the energy gain of the protons in each case. The accelerating electric field of the carbon case is higher than 
the gold case and the energy gain of the protons in the carbon case is higher too. It is similar to the theoretical lines 
(Fig. [ljc)(d)). Figure [Sfb) shows the energy spectrum of the generated protons for the gold case and the carbon 
case. The average proton energy, fave, is 38MeV in the gold case and 63MeV in the carbon case. The proton energy 
in the carbon case is 1.7 times higher than the carbon case. The proton energy spreads are 8.3% in the gold case and 
5.6% in the carbon case. Figures |3|c),(d) show the particle distribution and the x component of the electric field, E^, 
which accelerates the protons, around the target at t =50. The black solid lines show the E^ on the x axis. The value 
of Ex in the carbon case at the proton's position is 2.4TV/m and the gold case is 1.4TV/m. That is, the protons in 
the carbon case undergo a stronger acceleration. 

Figure [4] shows the average proton energy for some different materials of the first layer. The horizontal axis is for 
Zi/A, where Zi is the charge state of the ion and A is the atomic mass number. The vertical axis is for the proton 
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FIG. 5: (a) Initial shape of the target and the laser pulse with oblique incidence, (b) the interaction of the target and laser 
pulse, and (c) the exploding first layer and the accelerated protons (color scale), (d) The proton energy spectrum normalized 
by its maximum at t=100. (e) Distribution of carbon ions and protons (color scale), a 2D projection is shown looking along 
the z axis. 



energy which is normalized by the energy in the gold case. The higher energy protons can be obtained by using a 
larger ratio of Zi/A. That is, higher energy protons are obtained by using a high qi/rrii material in the first layer of 
a double-layer target. 

In Fig. |5]we show the case of an oblique incidence of the laser pulse onto the carbon target, showing the highest 
rate of proton acceleration in the above consideration. The laser pulse incidence angle is 30° [17] . The small hydrogen 
layer with the diameter of 4/3 A is shifted below the target center [16] (see Fig 5|a)). The other parameters are the 
same as described above. Figures |5|b),(c) show the electric field magnitude and distribution of ions. Half of the box 
is removed to reveal the internal structure of the electric field. We see that the carbon ions are distributed over a 
wider area due to the Coulomb explosion. Figure |5|d) shows the energy spectrum of the generated protons at t = 100. 
We obtain protons with an average energy, fave^ of lllMeV, an energy spread, A^^/fave? of 2.8%, and the number of 
generated protons of 2 x 10^. Figure [sje) is the 2D projection onto the {x^y) plane. 

In conclusion, proton acceleration driven by a laser pulse irradiating a double-layer target, consisting of some high-Z 
atom layer and a hydrogen layer, is investigated with the help of 3D PIC simulations. We find that higher energy 
protons are obtained by using a material with a high charge-to-mass ratio in the first layer of a double-layer target. 
As seen in our simulations, due to the strong Coulomb explosion occurring in such a material, the protons keep 
accelerating for a longer time. We show that three times higher energy protons are obtained, even using the same 
laser pulse, by using the optimum material for the first layer and selecting the optimum incidence angle of the laser 
pulse. 

The computation as performed using the PRIMERGY BX900 supercomputer at JAEA Tokai. 
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